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, Abstract. Two signatures of the quark-gluon plasma - strangeness 'enhance- 

ment' and J/tp 'suppression' - in nucleus-nucleus collisions are critically discussed. 
^ 1 A recently developed statistical coalescence model for J/tp production is pre- 

O ' sented. The measurements at the RHIC energies are crucial for disentangling the 

' different scenarios of J/tp formation. 

oo : 

' ■ 1. Strangeness "enhancement" ? 

r- : 

, The idea of strangeness enhancement as a quark-gluon plasma (QGP) signal in 

nucleus- nucleus (A+A) collisions was formulated a long time ago [1]. It was based on 
the estimate that the strangeness equilibration time in QGP is of the same order (« 10 
fm/c) as the expected life time of the fireball formated in A+A collisions. Thus in the 
case of QGP creation the strangeness is expected to approach its equilibrium value in 
QGP. This equilibrium value is significantly higher than the strangeness production 
in nucleon-nucleon (N+N) collisions. Strangeness production in secondary hadronic 
interactions was estimated to be negligible small (this appears to be not correct!). 
Qh, Therefore, if QGP is not formed, the strangeness yields would be expected to be 

much lower than those predicted by equilibrium QGP calculations. Thus at that 
time a simple and elegant signature of QGP creation appeared: a transition to QGP 
should be signaled by an increase of the strangeness production to the level of QGP 
equilibrium value. 

5^ ■ How to check this prediction? The actual study has been done in the following 

way. The strangeness to pion ratio, 

E, = W + <f +g> , (1) 
(vr) 

was measured and analyzed. Different ratios of this type, f s — (if ± )/(7r), 
(f2)/(7r), (K^) / (N p ) , ...,(ST)/N P , have been also studied (N p is the number of nucleon 
participants in A+A collision). One expected that all these ratios should increase 
strongly in A+A collisions if the QGP was formed. On the other hand, the strangeness 
to pion ratio (1) (as well as other specific ratios f s ) increases with collision energy y/s 
in N+N collision too. To reveal the specific strong increase of the strangeness/pion 
ratio in A+A collisions due to the QGP formation the strangeness enhancement factor 
was introduced: 

E^(y-s) 
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where Ef an d E^ N correspond respectively to A+A and N+N collisions at the 
same cm. energy y/s. Does the strangeness enhancement factor R s (2) become 
large if the QGP is formed? The confrontation of this expectation with the data 
was for the first time possible in 1988 when the results from S and Si beams at 
SPS and AGS were presented. The experiment NA35 reported that in central S+S 
collisions at 200 AGeV the strangeness to pion ratio is indeed 2 times higher than 
in N+N interactions at the same energy per nucleon. But even larger enhancement 
(R s is about of 3) was measured by E802 in Si+A collisions at AGS. Recent data on 
central Au+Au collisions at low AGS energies 4 + 10 AGeV completed the picture: 
strangeness enhancement is observed at all energies, and it is stronger at lower energies, 
i.e. the function R s (2) increases monotonously with decreasing of \/s. At the low 
AGS energies one does not expect a creation of the QGP and therefore a substantial 
strangeness enhancement is evidently of a different origin. Thus AGS measurements of 
strangeness enhancement larger than that at SPS show clearly that the simple concept 
of strangeness enhancement as a signal of QGP does not work. 

Let me return to the observable Ef A and its energy dependence shown in Fig. 1. 




Figure 1. Collision energy [F = (i/s — 2mjv) 3/ ' 4 /v / ^ 1 ' /4 ) dependence of 
strangeness to pion ratio (1) for central A+A collisions (closed points), N+N and 
p + p colisions (open points). The prediction of the statistical model of Ref. [2] 
is shown by solid line. A transition to the QGP is expected between the AGS 
(F ft 2 GeV 1 / 2 ) and the SPS (F ft 4 GeV 1 / 2 ) energies and leads to the non- 
monotonic dependence of the strangeness to pion ratio. At high collision energies 
the ratio saturates at the value characterestic for equilibrium QGP. 



The statistical model of the early stage [2] leads to the following predictions for 
strangeness production: 

a). A non-monotonic collision energy dependence of the strangeness to pion ratio (1). 
This is due to the fact that the phase transition is expected to occur at an energy 
where the strangeness to entropy (pion) ratio in the confined (hadron) matter is higher 
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than in the QGP. Therefore, a creation of the QGP in the energy region between the 
AGS and SPS would change an initial fast increase of the ratio E s (1) in the hadron 
gas by a decrease to the level expected in equilibrium QGP. 

b). Very similar strangeness to pion ratio is predicted for SPS, RHIC and LHC 
energies as strangeness/entropy ratio in the QGP is almost independent of temperature 
(collision energy). 

Preliminary SPS data in Pb+Pb at 40 and 80 AGeV and RHIC data in Au+Au 
at y/s = 200 GeV seem to be in agreement with the above conclusions (see Fig. 1). 



2. J/tp suppression and enhancement 

A standard picture of J/tp production in hadron and nuclear collisions assumes a 
two step process: the creation of cc pair in hard parton collisions at the very early 
stage of the reaction and the subsequent formation of a bound charmonium state. 
Matsui and Satz proposed [3] to use J/tp as a probe for deconfinement in the study 
of A+A collisions. They argued that in QGP the colour screening dissolves initially 
created J/tp mesons into c and c quarks which at hadronization form open charm 
hadrons. As the initial yield of J/tp is believed to have the same A-dependence as 
the Drell-Yan lepton pairs, the measurement of a weaker A-dependence of final J/tp 
yield (J/tp suppression) would signal charmonium absorption and therefore creation of 
QGP. The production of charmonium states J/tp and tp' have been measured in A+A 
collisions at CERN SPS over the last 15 years by the NA38 and NA50 Collaborations. 
There are two unambiguous consequences of the standard J/tp suppression picture 
[3]. First, the number of J/tp particles produced before the suppression should be 
directly proportional to the number of cc pairs, N cS . Second, when the energy y/s 
and/or the number of nucleon participants N p of the colliding nuclei increase, the J/tp 
suppression becomes stronger, i.e. for the measured J/tp yield, (J/tp), the ratio 

(3) 

iV CC 

decreases with increasing of ^fs and/or N p . 

Recently the thermal model [4] and the statistical coalescence model [5, 6] 
were suggested for the charmonium production in A+A collisions. This statistical 
description makes very different assumptions about the underlying physics which 
generates charmonium. In the statistical coalescence model [5, 6] the J/tp particles are 
produced by recombination of c and c at the hadronization stage, and the picture is 
much different from the standard suppression scenario. For large number of cc pairs, 
N c - >> 1, the multiplicity of J/tp due to recombination of c and c should be roughly 
proportional to: 

N 2 - 

m) ~ -f , (4) 

where V is the system volume. In this case the ratio R (3) is expected to increase 
with increasing of y/s and/or N p . We call this J/tp enhancement. 

The equilibrium hadron gas (HG) model describes the hadron yields measured in 
A+A collisions in terms of three parameters: volume V, temperature T and baryonic 
chemical potential [Ib- This model reproduces basic features of the data in the whole 
AGS-SPS-RHIC energy region describing successfully the hadron yields (see e.g. [7]). 
For the RHIC energies the temperature parameter T is expected to be similar to that 
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for the SPS energies: T = 170 ± 10 MeV. The baryonic chemical potential becomes 
small (hb < T) and decreases with the collision energy. 

The HG model assumes the following formula for the hadron thermal multiplicities 
in the grand canonical ensemble (g.c.e.): 



Nj = |J / ir<h> 




cxp ± 1 



(5) 



where nij, dj denote particle masses and degeneracy factors. The particle chemical 
potential fij in Eq.(5) is defined as 

fj,j = bjfi B + sjfis + Cjfi c , (6) 

where bj 7 Sj,Cj denote the baryonic number strangeness and charm of particle j. The 
baryonic chemical potential hb regulates the baryonic density of the HG system 
whereas strange /is and charm \ic chemical potentials should be found from the 
requirement of zero value for the total strangeness and charm in the system (in our 
consideration we neglect small effects of a non-zero electrical chemical potential). 

The total multiplicities N* ot in the HG model include the resonance decay 
contributions: 

N j ot = N i + E B (^^« • ( ? ) 

R 

where B(R — > j) are the corresponding decay branching ratios. The hadron yield 
ratios Nj 0t /N- ot in the g.c.e. are then the functions on T and \ib variables and are 
independent of the volume parameter V. 

For the thermal multiplicities of both open charm and charmonium states the 
Bose and Fermi effects are negligible, and nij » T. Therefore, Eq.(5) is simplified 
to: 



m,-T 



3/2 



Nj^djVe^^) exp^j. (8) 

The HG model gives the J /%p yield: 

NT/% = Nj^ + Ry/)Nv + R( X i)N Xl + R(X2)N X2 , (9) 

where Nj/^, N$i, N Xl , N X2 are calculated according to Eq.(8) and R{ip') = 0.54, 
R(Xi) — 0.27, R{X2) — 0.14 are the decay branching ratios of the excited charmonium 
states into J/ij). 

In the canonical ensemble (c.e.) formulation (i.e. when the requirement of 
zero "charm charge" of the HG is used in the exact form) the thermal charmonium 
multiplicities are still given by Eq.(8) as charmonium states have zero charm charge. 
The multiplicities (8) of open charm hadrons will however be multiplied by an 
additional 'canonical suppression' factor (see e.g. [8]). This suppression factor is 
the same for all individual open single charm states. Therefore, if No is the total 
g.c.e. multiplicity of all open charm and anticharm mesons and (anti)baryons, the c.e. 
value of the total open charm is equal to: 

«r = n MM , (10 , 

where Io,h are the modified Bessel functions. To find No we use Eq.(8) for thermal 
multiplicities of the open charm hadrons in the g.c.e. and take the summation over 
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all known particles and resonances with open charm [9]. The canonical suppression 
factor Ii(No) / Io(No) in Eq.(lO) is due to the exact charm conservation. Therefore, 
the baryonic number, strangeness and electric charge of the HG system are treated 
according to the g.c.e. but charm is considered in the c.e. formulation where the exact 
charge conservation is imposed. For No << 1 one has Ii(No) / lo(No) — No/2 and, 
therefore, the c.e. total open charm multiplicity is strongly suppressed in comparison 
to the g.c.e. result. At the SPS energies the c.e. suppression effects are important 
for the thermal open charm yield even in the most central Pb+Pb collisions. These 
suppression effects become crucial when the number of participants N p decreases. 
Note that for No << 1 the multiplicities of the open charm hadrons are proportional 
to V 2 in the c.e. HG formulation (instead of V in the g.c.e.). 

The statistical coalescence model (SCM) [5, 6] assumes that the charmonium 
states are formed at the hadronization stage. This is similar to the thermal model of 
Ref. [4] . The thermal model [4] predicts that the J/tp to tt ratio is independent of yfs 
and N p at high collision energies. This is because both (J/tp) and (w) are proportional 
to the system volume and they both depend only on the hadronization temperature, 
Tjj = 170 ± 10 MeV, which is expected to be independent of yfs and N p at high 
collision energies. However, in the SCM the charmonium states are produced via a 
coalescence of created earlier cc quarks at the early stage of A+A reaction by the hard 
parton collisions. One needs then an additional parameter j c [5] to adjust the thermal 
HG results to the required number of N c -. This is analogous to the introduction of the 
strangeness suppression factor -f s < 1 [10] in the HG model, if the total strangeness 
observed is smaller than its thermal equilibrium value. We find j c > 1 so that the 
open charm hadron yield is enhanced by a factor j c and charmonium yield by a factor 
7^ in comparison with the equilibrium HG predictions. The c.e. formulation of the 
SCM is [6]: 

= \lcNo J^J^ + ll «H . (ID 

where Nh is the total HG multiplicity of particles with hidden charm. Note that the 
second term in the right-hand side of Eq.(ll) gives only a tiny correction to the first 
term, i.e. most of the created cc pairs are transformed into the open charm hadrons. If 
N c - » 1 one finds from Eq.(ll) that -f c N Q » 1, therefore, Ii(j c N )/Io{lcN ) -> 1, 
i.e. the g.c.e. and c.e. results for the open charm coincide. In this case Eq.(ll) is 
simplified to [5]: N C c — "f c N Q /2 + j 2 ^h- This happens for central Au+Au collisions 
at upper RHIC energy yfs = 200 GeV. For lower RHIC energy yfs — 56 GeV the value 
of N C c could be close to (or even smaller than) unity so that the c.e. suppression effects 
for the open charm are still important. Note that for the non-central A+A collisions 
the c.e. suppression effects become evidently stronger, hence their consideration are 
necessary to study the N p dependence of the charmonium production even at upper 
RHIC energy 

Eq.(ll) will be used to find the charm enhancement factor j c and calculate then 
the J/tp multiplicity: 

(J/V) = 7c 2 N!j% , (12) 

where Nf^ is given by Eq.(9). Note that T = 170 180 MeV leads to the value of 
the thermal ratio, 

Ml= (j*tL) m ^(- ^- m '/* ) ^0.04^0.05,(13) 



(J/tp) \m 



Strangeness and charm signatures in A+A collisions 



6 



in agreement with data [11] in Pb+Pb collisions at SPS for N p > 100. This fact 
was first noticed in Ref. [12]. Recent results [13] on transverse mass spectra of J/tp 
and tp' mesons in central Pb+Pb collisions at 158 AGeV also support a hypothesis of 
statistical production of charmonia at hadronization (see Ref. [14]). 

The number of directly produced cc pairs in the left-hand side of Eq.(ll) should 
be estimated in the pQCD approach and used then as the input for the SCM. The 
pQCD calculations for cc production cross sections were first done in Ref. [15]. For 
the cross section o~(pp — > cc) of the charm production in p+p collisions we use the 
results presented in Ref. [16]. This leads to the value of cr(pp — > cc) = 0.35 mb at 
y/s = 200 GcV and the -^/i-dependence of the cross section for y/s = 10-7- 200 GeV is 
parameterized as: 

a( PP ^^=a .(l-^Y (14) 



M 

with cr = 3.392 fjb, M = 2.984 GeV, a = 8.185 and = 1.132. 

The number of produced cc pairs in A+A collisions is proportional to the number 
of primary N+N collisions, N££, which in turn is proportional to Np^ 3 [17]: 

N c - C = N^(N p ) = C a(pp cc) , (15) 

a NN 

where cr™f} = 30 mb is the inelastic N+N cross sections, C = 11 barn -1 . 

The results of the SCM can be studied analytically in the limiting cases of small 
and large numbers of N c -. For N cZ « 1 one finds from Eq. (11) that 7^ ss 4N C c/Nq 
and then from Eq. (9) 

2^*0* 4/V to * 1 1 

This behavior is similar to the standard picture of the J/tp suppression: the ratio 
R decreases with increasing both y/s and N p . Therefore, the SCM predicts the J/tp 
suppression at the SPS energy. This energy is still too "low" as N C c < 1 even in the 
most central Pb+Pb collisions at 158 AGeV. This behavior is however dramatically 
changed at the RHIC energies [18]. In most central Au+Au collisions (N p = 2A) 
at y/s = 200 GcV the expected number of N C c is essentially larger than unity. For 
N c - » 1 one finds from Eq. (11) that 7 C w 2N cZ /N and then from Eq. (9) 

ll N *J% _ 2N C - N c - N c - /TN/3-1/2 w l/3 / 17 x 

where (3 = 1.132. Eqs. (16,17) reveal a remarkable prediction of the SCM: the J/tp 
to N c - ratio decreases with both ^fs and N p (J/tp suppression (16)) when N c ^ << 1 
and it increases with both y/s and N p (J/tp enhancement (17)) when N c ^ » 1. 
The measurements in Au+Au collisions at RHIC give a unique possibility to check 
simultaneously both these predictions: at y/s — 56 GeV and N p = 100 the expected 
value of N cZ from Eq. (15) is N c - = 0.2 << 1, but at = 200 GeV and N p = 2A 
one expects from Eq. (15) that N C c = 10 >> 1. Therefore, changing of y/s and N p 
in Au+Au at RHIC one could observe both the J/tp suppression and enhancement 
behaviors. These limiting behaviors are smoothly connected in the intermediate region 
of and N p where N c - == 1. The results of the SCM for the J/tp to N C c ratio are 
presented in Figs. 1 and 2. Both the suppression (the dashed line in Fig. 3) and 
enhancement (the solid lines in Figs. 2 and 3) behaviors are clearly seen. 
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Figure 2. The energy dependence of the J/ip to iV c g ratio in central (N p = 2A) 
Au+Au collisions. Points are the predictions of the SCM for the RHIC energies: 
y/s = 56,130,200 GeV. The ratio R increases by a factor of about 3 (J/tp 
enhancement) in the region of RHIC energies yfs = 56 4- 200 GeV. 




Figure 3. The iVp-dependcnce of the J/tp to N c - ratio. The lines are the 
predictions of the SCM. The dashed line corresponds to y/s = 56 GeV and shows 
the J/tp suppression behavior. The solid line corresponds to yfs = 200 GeV and 
shows the J/tp enhancement. 
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3. Conclusions 

• The dcconfincmcnt phase transition is expected to occur at the collision energies 
between AGS and SPS where the strangeness to entropy (pion) ratio reveals a 
non-monotonic (or kinky) dependence on the collision energy [2] (see Fig. 1). 

• Statistical hadronization of the QGP is probably an important source of J/tp 
production [4]. This fact would open a new look at J/tp 'suppression' signal of 
the QGP. The statistical coalescence model [5, 6] of the J/tp production predicts 
that the J/tp suppression in peripheral Au+Au collisions at lower RHIC energy 
should be changed into the J/tp enhancement [18] in central Au+Au collisions at 
the upper RHIC energy (see Figs. 2 and 3). 
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